First passage time for subdiffusion: The nonadditive entropy approach 

versus the fractional model 
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We study the similarities and differences between different models concerning subdiffusion. More 
particularly, we calculate first passage time (FPT) distributions for subdiffusion, derived from 
Greens' functions of nonlinear equations obtained from Sharma-Mittal's, Tsallis's and Gauss's non- 
additive entropies. Then we compare these with FPT distributions calculated from a fractional 
model using a subdiffusion equation with a fractional time derivative. All of Greens' functions 
give us exactly the same standard relation ((Ax) 2 ') = 2D a t a which characterizes subdiffusion 
(0 < a < 1), but generally FPT's are not equivalent to one another. We will show here that the 
FPT distribution for the fractional model is asymptotically equal to the Sharma-Mittal model over 
the long time limit only if in the latter case one of the three parameters describing Sharma-Mittal 
entropy r depends on a, and satisfies the specific equation derived in this paper, whereas the other 
two models mentioned above give different FTPs with the fractional model. Greens' functions ob- 
tained from the Sharma-Mittal and fractional models - for r obtained from this particular equation 
- are very similar to each other. We will also discuss the interpretation of subdiffusion models based 
on nonadditive entropies and the possibilities of experimental measurement of subdiffusion models 
parameters. 

PACS numbers: 02.50.Ey, 05.40.-a, 66.10.C-, 05.70.Ln 



I. INTRODUCTION 

Over about the last 20 years the anomalous diffusion 
process has been observed in many physical systems. Si- 
multaneously, various theoretical models of this process 
have been put forward (see [l|-|3[ and the references cited 
therein). It is worth considering what the issue of anoma- 
lous diffusion is and what definition of this process can 
be taken into account. We note that the situation is 
different here compared to the normal diffusion process. 
Namely, in the latter case the different models produce 
results which are equivalent to each other. For exam- 
ple, the stochastic approach provides the same normal 
diffusion equation as entropy formalism, and their fun- 
damental solutions appear to be Gaussian functions with 
their second moment proportional to time 



((Ax) 2 ) = 2Dt 



(1) 



Consequently, the question: 'What is the definition of 
normal diffusion?' can be answered in a few equiva- 
lent ways, such as with the process described by the 
Langevin equation with white noise, the standard Wiener 
process, the random walk of a particle where the prob- 
ability distributions of its step length and the waiting 
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time to take its next step have finite moments, or the 
process described by a probability density which maxi- 
mizes the Boltzmann-Shannon entropy. Although from a 
mathematical point of view these definitions are not ex- 
actly equivalent to each other, physicists usually treat 
these definitions equivalently. Anomalous diffusion is 
a process which qualitatively differs from normal diffu- 
sion, so all of the above mentioned definitions are not 
fulfilled. However, such treatment causes ambiguity in 
the definitons. There arises a problem with the defini- 
tion of anomalous diffusion as a process which denies 
normal diffusion definition. The Wiener process is re- 
placed by the frational one (consequently the Langevin 
equation is changed), whithin the Continuous Time Ran- 
dom Walk (CTRW) formalism at least one of the proba- 
bility distributions describing a single particle jump has 
infinite moments, and nonadditive entropies are used in- 
stead of Boltzmann-Shannon entropy. The anomalous 
diffusion is frequently defined by its interpretation as a 
non-Markovian random walk which is described by non- 
Gaussian probability distribution. However, this defini- 
tion is limited to the stochastic processes and its relation 
to thermodynamics or deterministic models is not obvi- 
ous. To find a more general definition one sholud assume 
that anomalous diffusion is characterized by a special pa- 
rameter (in the following denoted by a) which is a 'mea- 
sure' of how far the anomalous diffusion process is from 
the normal diffusion one. Since within the CTRW for- 
malism the anomalous diffusion model provides us with 
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the probability distributions with the second moment to 
be nonlinear of time 



((Ax) 2 ) = 2D a t a 



(2) 



where a > 0, a^l, the statement that the model which 
provide the relation ([2]) can be treated as an anomalous 
diffusion model has been rather widely used. The param- 
eter a plays a different role in the models; it is related to 
the fractional derivative order in the anomalous diffusion 
equations or it controls a measure of nonadditivity of the 
entropy. 

The most used anomalous diffusion models seems to be 
the CTRW model and the models derived from nonaddi- 
tive entropies. Throughout this paper we will refer to the 
model based on CTRW as a fractional model. CTRW 
provides the linear anomalous diffusion equation with 
the fractional-order derivatives 2-4], while the models 
based on the nonadditive entropies give nonlinear differ- 
ential (or integral-differential) equations with derivatives 
of a natural order [5l4l4|. The simplest stochastic inter- 
pretation of the anomalous diffusion seems to be found 
within CTRW models, where the random walker waits an 
anomalously long time to make its next step, for exam- 
ple in the transport process of sugars in gel (for subdiffu- 
sion) or where the step length can be anomalously long 
with a relatively high probability, for example a random 
walk in a turbulent medium (for superdiffusion). For the 
nonadditive entropy model, a physical meaning of the 
anomalous diffusion equation is manifested mainly in its 
stationary version, namely, the stationary solution of the 
equation maximizes nonadditive entropy. However, as is 
shown in (lEj . there is a non-Markovian process which 
provides the relation (fTJ). Thus, the stochastic definition 
of anomalous diffusion is, in some situations, in contra- 
diction with the ones based on the relation @. 

We note that in a system where the relation ([2]) is valid, 
other functions, which can be measured experimentally, 
of a type f(t) ~ t a / 2 describe the system, such as the 
time evolution of near membrane layer thickness flU \vf\ 
and the time evolution of the reaction front in the subdif- 
fusive system with chemical reactions [l8l . [l9j . Let us also 
note that there are models which do not have fully sat- 
isfactory stochastic interpretations yet, as, for example 
the anomalous diffusion process described within nonad- 
ditive entropy formalism; but in these models we can also 
find the power-like important characteristic of the sys- 
tem identical with the ones found within stochastic mod- 
els. The stochastic models are based on the assumptions 
which simplify the problem which the experimental data 
is not capable of confirming, as the assumptions con- 
cerning random walk of a particle for a separable case 
of CTRW formalism. On the other hand, the stochas- 
tic interpretation of various diffusion models is needed. 
The attempts to find stochastic interpretation of nonad- 
ditive entropy formalism have mainly been made using 
the modified Langevin equation for the description of the 
anomalous diffusion process which is simultaneously de- 
scribed by a nonlinear differential equation. The result is 



that a random force occurring in the Langevin equation 

tends on the solution to the nonlinear equation [20l — 
This situation can be interpreted as the existence of 
feedback between a system and a random force, whose 
interpretation is — at least in our opinion — not obvi- 
ous. We remark here that such nonlinear equations were 
also derived from the master equation [25| . Recently, 
we found a new stochastic interpretation of sub-diffusion 
as a 'long memory diffusion' described by the genera- 
zlied linear Langevin equation in a system with external 
Gamma-type noise [26]. The anomalous diffusion equa- 
tion derived within this new model perfectly coincides 
with the one derived within the nonadditive Sharma- 
Mittal entropy formalism. 

The above considerations lead us to take the following 
definition of the anomalous diffusion: An equation whose 
fundamental solution (Green's junction) G(x,t;xo) pro- 
vides the relation (0) where 



<(Ax) 2 ) 



(x — xq) 2 G(x, t; xo)dx 



can be treated as an equation describing the anomalous 
diffusion process. Throughout this paper we have as- 
sumed that we are dealing with one-dimensional homo- 
geneous systems without any external fields and convec- 
tive flows; here D a is the anomalous diffusion coefficient 
measured in the units m 2 / s a and a is the anomalous dif- 
fusion parameter; < a < 1 for sub-diffusion, a > f for 
superdiffusion. Thus, anomalous diffusion is controlled 
by two parameters: D a and a. We would like to add 
that the anomalous diffusion coefficient is sometimes de- 
fined differently, for example the subdiffusive coefficient 
D a is often defined as 0,0] D a — T(l + a)D a . Although 
Greens' functions derived from the equation mentioned 
above provide Eq. ([2]), other important characteristics, 
such as a first passage time distribution, are different. 

First passage time (FPT), which is one of the most 
important characteristics in normal and anomalous diffu- 
sions [27l - l29j , was studied for anomalous diffusion mainly 
within the fractional model 2, 3.j30l-l4l| as well as in the 
lattice and fractal medium [42|, |43| . The FPT is defined 
as the time that the random walker takes to reach a tar- 
get located in xm for the first time, from the starting 
point xo. The FPT has been used to describe real physi- 
cal processes such as animals searching for food 44j , the 
passage of DNA molecules through a membrane channel 
[30j | and the spreading of disease [HI etc. Moreover, the 
FPT distribution can be used to calculate other char- 
acteristic functions which are measured experimentally, 
such as the time evolution of an amount of a substance 
leaving a sample or the substance flux flowing through a 
sample surface. 

In this paper we study subdiffusive systems. We de- 
rive the FPT distributions for subdiffusion equations de- 
rived form Sharma-Mittal, Tsallis and Gauss nonaddi- 
tive entropies and compare them with the result obtained 
from the fractional model. We sholud add here that we 
will pay particular attention to the nonlinear differential 
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anomalous diffusion equations derived from nonadditivc 
entropies and we will not discuss the interpretation of the 
entropies. Our strategy is as follows: we adapt Greens' 
functions presented in their generalized forms in Frank's 
book [5[ into special forms so that each of them exactly 
satisfies the relation ([2]). Next, we calculate the FPT 
distribution functions and compare them with the one 
obtained from the fractional model. We derive consis- 
tency conditions in the functions calculated from various 
models over the long time limit. We briefly discuss the 
similarity of Greens' functions and the interpretation of 
the parameters occurring in nonadditive entropies and 
their connection with the parameters D a and a, which 
are measured experimentally. This method would allow 
us to measure nonadditivc entropy parameters. 



combined with the continuity equation dC/dt = —dJ/dx 
For the normal diffusion case L(P) ~ 1/-P- There is 
no obvious choice of the function L(P) for anomalous dif- 
fusion. The anomalous diffusion equation derived within 
nonadditivc entropy formalism reads [5] 

^^.[p,™, (4) 

where 1 and r>l/3,r = qfor the Tsallis case and 
r = 1 for the Gauss case, the index i denotes the symbol 
identifying entropy, Qi denotes the fluctuation strength, 
and 




*t[P] = 1 , 



II. THE METHOD 

Here we will present the functions and equations which 
define Green's function and distribution of first passage 
time used in our considerations. The description of non- 
linear anomalous diffusion equations derived from non- 
additive entropies and their analytical solutions (Greens' 
functions) are based on Frank's book [5] . 



A. Anomalous diffusion equations 

1. Nonadditive entropy formalism 
The Sharma-Mittal entropy is defined as 

,(?-l)/(r-l) 



oSM 
q,r 



[P] = 



l-(JP r dx) 



(3) 



where q, r > 0, q, r 7^ 1, P denotes a probability of find- 
ing a particle at the point x at time t. Gauss entropy 
obtained in the limit r — ¥ 1~ reads 



f 



S^[P] 



where q > 0, q ^ 1. The Tsallis entropy can be obtained 
from ([3]) putting q = r. For two statistically independent 
systems A and B the Sharma-Mittal entropy satisfies the 
following equation 



S™(A + B) = S»™(A) + 



SMi 



S s .t\B) 

<7 r (A)Sg*? (B) 



q. r 



For q 7^ 1 one deals with nonadditive entropy [46l. l47j . 

There are two main ways to derive anomalous diffu- 
sion equation from nonadditive entropy. Within the first 
one the stationary state is generated by means of a max- 
imum entropy condition under conditions which assume 
that the g-moments are finite. Within the second method 
the flux J defined by the equation J = L(P)(6S/SP) is 



2. Fractional model 

Within the separable case of CTRW formalism it is 
assumed that a particle takes its next step of a length 
p after time r, where both are independent random 
variables. For subdiffusion the probability distribution 
uj(t) is of 'thick tail', w(r) w — r Q /i 1+Q T(— a) for a suf- 
ficienty large time (the mean value of u> equals infin- 
ity), whereas X(p) is the Gaussian distribution, A(p) = 
exp (— p 2 /2cr 2 ) /V27TCT 2 . Under the above assumptions 
one obtains the linear subdiffusion differential equation 
with the Riemann-Liouville fractional derivative 0tJ] 



dP(x,t) _ - & 



d 2 P{x,t) 



dt 



dt a dx 2 



(5) 



where D a — a 2 /r a , the Riemann-Liouville fractional 
derivadivc is defined as follows for a > 



d a f(t) 



1 



(t-f) 



f\n—a — l 



f(t')dt' 



dt a T(n - a) dt n J 
where n is a natural number fulfilled a < n < a + 1 . 

B. Green's function 

Green's function (GF) is defined here as a solution to 
the appropirate diffusion equation with the initial condi- 
tion 

G(a;,0;a;o) = 5{x - xo) , 

5 denotes the delta-Dirac function. When particles are 
independently transported and all of them start their 
movement at xq at the initial moment t = 0, Green's 
function can be interpreted as a concentration profile of 
the particles normalized to one (i.e. divided by the num- 
ber of particles). Another interpretation of the GF is 
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that it is treated as a probability density of finding a 
particle at point x at time t under the condition that its 
initial position is xq. We should add here, that there is 
a stochastic interpretation of subdiffusive movement of 
a particle described by a nonlinear differential equation, 
which assumes that the particle is transported indepen- 
dently of other particles [26| . For the unrestricted sys- 
tem the GF satisfies natural boundary conditions which 
require the disappearance of the function at an infinite 
distance from the initial position G(x,t;xo) — > 0, when 
x — » ±00. 

Green's function for the nonlinear diffusion equation 
(m >0,m^l) 

dP(x,t) d 2 P m {x,t) 



01 



dx 2 



known as the Barenblatt solution, reads [48 



G(x, t; xq) = t 



1 - 



k(m — l)\x — xq\' 
2mt 2k 



(6) 

where k = l/(m + 1), {u} + = max{w, 0}. 

In our paper we use Green's function for the system 
with a fully absorbing wall. The commonly used bound- 
ary condition at an absorbing wall reads 

GabeixAU^Xo) = . 

Due to the interpretation of GF and the symmetry argu- 
ments, the GF for the normal and subdiffusive systems 
with a fully absorbing wall can be found through the 
means of the method of images, which for x,xq < xm 
gives 

G abs (a;, t; x ) = G(x, i; x ) - G(x, t; 2x M - x ) . (7) 



1. Nonadditive entropy formalism 

Greens' functions of Q for the Sharma-Mittal and 
Tsallis models take the form of Eq. ([6]) and read [5| 



GsM(x,t;x ) 
= D SM (t) 



CsM(t) , ,w x 2 

1 (r- l)(x-x ) 



(8) 



where r > 1/3, r, q^ 1, q > 0, 
GT(x,t;x{)) = 



D T {t) 



{q-l){x-x Y 



where q > 1/3, q ^ 1, 
G G (x,t;x Q ) = D G (i)exp 



C G (t) 



(x - XqY 



(9) 



(10) 



The functions occurring in (l8|- (lT0l) are defined as 
D SM (t) ' 



2r(l + q)QsMK Tt q\z r \ t 



D T (t) = 



D G (t) 



2q{l + q)Q T \z q \U 



3 (9-l)/2 



2tt(1 + q)Q G t 



(ID 



(12) 



(13) 



where 



r(r/(r-l)) 



f=Tr((3r-l)/(2(r-l))) ' T > 1 > 

/7F , r = 1 , 

~r((l+r)/2(l-r)) 
1-r r(l/(l-r)) 



(14) 



1/3 < r < 1 . 



and 



r=l, 



CsmW =2(z r % f (i)) 2 , 
C r (t) = 2(z q D T (t)) 2 , 
C G (t) = 2n{D G {t)) 2 . 



(15) 



(16) 



The second moment of Greens' functions presented above 
can be calculated according to the formulae 

2 1 



(AxY(t) 



sm 3r -1Csm(*)' 



(17) 



:(*) 



i8. Fractional model 



For Eq. ([5]) the Green's function is (in the following, 
the index F corresponds to the fractional model) 



n i . s 1 , I . k-a;o| 
GF(x,t;xoJ = t=/q/2-i.q/2 t; 

2VA^ V , 

where 

1^ 1 / a\ k 

UAt; a) - 2^ T( _ k0 _ m 



k=0 



(18) 



(19) 



/3, a > 0, the function / can also be expressed in terms 
of the Fox function 49]. We should add here, that the 
methods of solving Eq. (|S|) are presented, among others, 

in @, m, 
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C. First passage time 

Let us assume that a particle is located at Xq at the 
initial moment t = 0. The time when the particle reaches 
the point xm for the first time is a random variable de- 
scribed by a probability density of F, calculated accord- 
ing to the formula for t > 



F(t;x ,XM) 



dR(t;x ,x M ) 
dt 



(20) 



for t < we put F(t; xq,xm) = 0, and where R(t; xq, xm) 
denotes the probability of finding the particle at time t 
starting from xq in the system with a fully absorbing wall 
located at xm (in the following we assume that xq < xm) 



R(t;x ,x M ) 



G ahs (x,t;x )dx . (21) 



The cumulative function of F is 1 — R(t; xq,xm)- The 
mean first time (T) is defined as 

/>oo 

(T) = / tF(t;x ,XM)dt . 
Jo 

Assuming that tR(t;xo,X]\i) — » when t — » oo, one ob- 
tains 



sets using the FPT distributions. Our further consider- 
ations are based on the following assumptions: the first, 
since the models describe the same subdiffusion process, 
all of Greens' functions should provide the same relation 
which defines subdiffusion ([2]); the second, the parame- 
ters a and D a are measured experimentally. The exam- 
ples of such measurements, where the empirical results 
were compared with theoretical functions derived within 
the fractional model, are presented in [lfl [l7|> HH ■ Taking 
into account Eqs. ©, (HI]), ([Tjl]), ([13]), QI]), and CE7]) we 
come to the conclusion that relation @ will be satisfied 
by all of Greens' functions Gi(x, t; xq) when 



q = --l 

a 



(23) 



for all models, while the fluctuation strength must be cho- 
sen for each model separately. The fluctuation strength 
Qi, which in general depends on q and r, plays the key 
role in expressing the parameters of the models based on 
nonadditive entopy by a and D a . When q = r = 1 we 
are dealing with normal diffusion and then Qi is identi- 
fied as the normal diffusion coefficient, D\ — Qi for each 
model. Taking into account (|2"3"]l they read 



lV« 



!SM 



a[2D a (3r- X 
^K r , 2la _ x \z r \^l^ 



(24) 



(T) = / R(t;x ,x M )dt 
Jo 



(22) 



l/a 



a 2 [2D a (6/a - 4)] 
4(2 -a)!^!^ 1 - 1 /") 



(25) 



III. GREENS' FUNCTIONS GENERATED 
RELATION © AND FIRST PASSAGE TIME 
DISTRIBUTIONS 

Greens' functions (]8]), (|9]) and (fTOj) . presented in the 
previous section, depend on the parameters q (which can 
be interpreted as a measure of entropy nonadditivity) 
and the fluctuation strength Qi, here the index i de- 
notes the model, i = SM, T, G for the Sharma-Mittal, 
Tsallis and Gauss models, respectively. Sharma-Mittal 
Green's function additionally depends on the parameter 
r. Green's function for the Tsallis and Gauss models can 
be treated as specific cases of the Sharma-Mittal one, 
using q — r for the Tsallis model and having as a limit 
r — > 1~ for the Gauss model. However, retaining the 
commonly used terminology, we consider these functions 
separately. 

We look at two sets of subdiffusive models having dif- 
ferent physical origins. The first contains the models 
derived from nonadditive entropy and the second con- 
tains the fractional model. In general, both of these sets 
can describe the same physical processes. However, pro- 
cesses also occur which can be described by models from 
one set alone, the other set of models cannot be appli- 
cable in describing such a process (this problem will be 
briefly discussesd in the Final Remarks) . We are going to 
find the accordance conditions between models from both 



a(2D c 



\l/a 



2(l-l/a) 



(26) 



In the following subsections, we apply Eqs. (|2"3"]l - (f2"6"|) 
to eliminate the parameters q and Qi from Greens' func- 
tions ([8]), ([9]) and (|T0]) and then we use Greens' modi- 
fied functions to calculate the functions Fi(t; Xo,Xm) and 
Ri(t;xo,XM) from (12"0]) and (f2"Tj) . In the last subsection 
we will find these functions for the fractional model. 



A. The Sharma-Mittal model 

Using Eqs. ©, (|2"3"]). and dM]), Green's function for 
the homogeneous unrestricted system, provided by the 
Sharma-Mittal entropy model, reads 



GsM{x,t:x ) = 
1 



^2D a (3r - l)t a \z r \ 



(r - l)(x - Xq) 2 

2D a (3r - l)t a 



(27) 



where r > 1/3 and r / 1. The function has different 
properties depending on the value of the parameter r. In 
the following, we consider the cases r > 1 and 1/3 < r < 
1 separately. 
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FIG. 1. Green's function for SM model in the case of r > 1 
with finite support, W = 154.7 for r = 1.125, 79.3 for r = 2, 
66.4 for r = 5, 63.7 for r = 10, 61.8 for r = 50 and 61.7 for 
r = 100, here xq = 0. Values of the parameter r are given in 
the legend. 



1. The case of r > 1 

The function (f2"T|) has a finite support for r > 1 , so the 
probability of finding the particle differs from zero only 
in the interval 

xe (x -W(t),x + W(t)) , 



where W(t) = Bt a / 2 , B = 



_ /2D a (3r-l) 



(see Fig. H). The 



boundaries of the interval move with the speed v g given 
by the relation 



dW(t) 
dt 



aB 



2tl-«/2 ■ 



(28) 



The finiteness of v g ensures that there is a minimum 
time T X2tXl of the passing of the particle from the point 
x\ to a?2, which is given by the relation 



Tr 



(x 2 - Xi) 2 
B 2 



l/c 



(29) 



Using 0, d2TJ), ([211) and ((29]), we obtain 

RsALr>l(t',X ,X M ) = 

= e(T XMtX0 -t) + e(t-T XMiX ^ 

a/2 



\z r \y/r ~ 1 



t 



2F1 



1 -1 3 /T s 



2' r- 1 ! 2' 



(30) 

where is the Heaviside function and 2-P1 [a, c; 2] de- 
notes the hypergeometric function, 



2 Fi [a, 6; c; z] 



\ - (a)„(b)„ z n 

h (<=)» nl 



(31) 




FIG. 2. Green's function Gsm for 1/3 < r < 1. Values of the 
parameter r are given in the legend. 



(a)„ is the Pochhammer symbol, (a) n = T(a + n)/T(a). 
From (EB and CH we obtain 



FsM,r>i(t; X , Xm) 



(32) 



|z r |\A- - 1 



1 - 






V * / J 



41+01/2 

2. ITie case o/l/3 < r < 1 



For this case Green's function is unrestricted. The tails 
for a specific given t reads G(x,t;x ) ~ l/x 2 ^ 1_r ^ when 
a; — » ±00. The graphs of Greens' functions for this case 
are presented in Fig. [5J By repeating the procedure 
presented above we obtain 



RsM,r<l(t;Xo,XAl) 



a/1 - r\z r \ 



(XM - Xq ) 2 



2^ Q (ff^) + (^Af-a:o) 2 



1 1 - 3r 3 



(x M - ^o) 2 



2' 2(1 - r) 2 2A^(^i) + (x M ~ x )< 



and 

FsM,r<l(t', Xq, X M ) 
_ a\x M - X \ 



(33) 
(34) 



2D 



( 

°< \ 1-r J 



c(l + r) 
t2(l-r) ' 



|z r |Vl - r 



2D 



a (t5f) t" + ( X M - Xq) 2 



Eqs. (ISU)) . (I3lj) and ([33]) provide the same formula 
for long times, namely for t 3> tsM, where t$M = 



7 



(X M -X ) 2 



l/a 



2D Q |(3r-l)/(l-r)| 

RsM(t;x ,x M ) = 
for r > 1/3 and r / 1. 



we obtain 

V2\x m -xq\ 1 
= ^D a (3r-l)\z r \ t°/ 2 ' 



(35) 



For long times t^$> ty, where 



(2-aQ(x M -z ) 2 
4D Q (3-2a) 



Eq. (|38p can be approximated as follows 

\x M - ^01^2(1 -a)T (2(1^7) 



l/a 



R T (t;x ,x M ) 



fa /2 



B. The Tsallis model 

As we mentioned previously, Tsallis Green's function 
can be treated as a specific case of the Sharma-Mittal 
one. The results presented in the previous section are 
valid for the Tsallis model in which 



r = 1 

a 



(36) 



Let us note that the Tsallis model corresponds to the 
case of r > 1 for subdiffusion. As a formality, we present 
the functions derived in Sec. IIII A 11 taking into account 
Eq. 

G T (x 7 t;xa) = 



^2D a (3-2a)\z 2/a „ 1 \t<*/ 2 



where 



z 2/a-l 



(1 - a)(x - XqY 
(6 - 4:a)D a t a 



(271=^) 



(37) 



on 



R T (t;x ,x M ) 



&{T XMiX0 -t)+O(t-T XMiX0 ) 



\2(l- a )) V 



2r ; 



2a; » 



> / 2 



2F1 



-a 3 I T x 



2' 2(l-o)' 2' 



(38) 



where 



(2 - a)(afAf - x Q ) 2 



4D a (3- 2a) 
and the distribution of the FPT is 

F T (t;x ,x M ) = 

r (5TJ L 7iy) 



l/a 



(39) 



&(t-T XMtX0 )- 



^(277=^ 



fl+a/2 



C. The Gauss model 

Green's function for the Gauss entropy model can be 
obtained from (1271) in the limit r —> 1~ and reads 



G G (x,t;x ) = 



1 



exp 



(x - X ) 2 
4D n t a 



( x M - a^o 



2y/wD a t a 
From (J3TJ) and (gDJ) we obtain 

^ G (^o,XM) = erf vy __ 

where erf(u) is the error function 



t , S 2 f -u 2 , 2 ~ (-1)« M 2 "+! 

erf(u) = — = / e u du = — = > ^— - — — — 
V^io v^fr/ (2n+l)n! 

Using ((201) and 6D we g et 



F G (i;z ,x M ) = 1 M ^-°' ^-^e~ 



2^D~ a t l + a ' 2 



(40) 



(41) 



(42) 



(43) 



For t ^S> to, where to — 
(|4*21 we obtain 

R G (t;x ,x M ) 



AD a 



l/a 



from (|4"Tj) and 



D. The fractional model 



The Green's function (fT5|) provides the relation ([2]) if 



D a = T(l + «)£>« . 
Applying Eqs. Q, ((IS) and ([21]) we get 



i? F (i; x ,x M ) = 1 - /_i ia /2 I i 
and from (1451) and (l20l) we obtain 



F F (t;x a ,x M ) = /o,a/2 



For £ 3> tp, where if = 
Eqs. ((Til), (HSJ), and (gS) give 



3-M - Xp 



XM - Xp \ 

J 

|a?M-a:Q|r(l-a/2) 



R F (t;X ,X M ) 



y/D a T(l + a)T(l - a/2) t a ' 2 



(44) 
(45) 

(46) 

1 2/a 

1 

(47) 
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FIG. 3. Comparison between F (t; xo, xm) for the Sharma- 
Mittal model for different r values and the other models. 




le+05 le+06 



log 10 (0 

FIG. 4. The same situation as in Fig. [3] but in the log-log 
scale for long times. 



We should add here, that the equivalent results to 
Eqs. (|4"5l) . (PfBl) and (H7)) were previously obtained by 
Barkai |4l| . 



IV. COMPARISION OF THE MODELS 

In this section we compare the functions obtained from 
the different models. Putting the functions (j30]l . (J33J) , 
dMJ), flU]), (05]) into (gU respectively, we get (T) = oo 
for all models. 

The probability densities of FPT are presented in Figs. 
[3] and HJ In both graphs the functions obtained for the 
fractional model are represented by a dashed line, for the 
Tsallis model - by a solid one with empty circles and for 
the Gauss model - by a dotted line; the other lines are 



FIG. 5. Comparision of the numerical solutions of Eq. (|49p 
(circles) and the plot of the function (|50|l (solid line). 



assigned to Sharma-Mittal functions with various r. All 
graphs are prepared for xq — and xu = 100. The val- 
ues of the rest of the parameters are given in each figure 
separately (all quantities are given in the arbitrary chosen 
units). In Fig. [3] we have presented graphs for relatively 
short times calculated for all models studied in this pa- 
per, and in Fig. [4] we have presented these functions for 
long times (in the log-log scale). We can observe that for 
a very short t the function F$m is larger for a smaller r 
value (for r > 1 this function equals zero for t < T XMtXo ). 
For intermediate length times we observe the opposite 
situation, but the Gauss function is the largest now. For 
very long times, the situation again changes to the op- 
posite and the largest value takes the function for the 
smallest r value, whereas the Gauss function takes me- 
dial values. In the log-log scale (Fig. S|) we observe that 
for sufficiently long times, the tails of F(t; xq, xm) for all 
models and different r values are parallel, and accord- 
ing to the functions ([32]), ((34j) , <j39j), (gSJ), and (06]) there 
is F ~ \/t 1+a / 2 . Let us note that the functions become 
parallel for a different time, and this time depends on the 
values of the 'reference time' ti, i = SM,G,T, F. The 
tsM depends heavily upon the parameter r and for the 
cases presented in the graphs there are tsM — 192417 for 
r = 0.375, 19307 for r = 0.5, 7650 for r = 0.59, 288 for 
r = 1.125 and 3981 for r = 100, the other reference times 
are: t T = 1767, t G = 7173, and t F = 794. We would like 
to add that the graphs of the functions F$m for various 
1.125 < r < 100 are located between the functions corre- 
sponding to r = 1.125 and r — 100 (for graphical clarity 
these functions are not presented in graphic figures), for 
r > 100 all functions are very similar, and in practice 
they are difficult to distinguish. Let us note that some 
of the functions F calculated from the Sharma-Mittal 
model are very similar to the ones obtained from the 
fractional equation for sufficiently long times for some r 
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- r=0.64, D a =5 
0,008 - g=0.3, t=10 



0,002 - 




FIG. 6. Comparison between Greens' functions for the 
Sharma-Mittal (solid lines) and the fractional (dashed lines) 
models for different times given in the legend for r given by 
ED. 



FIG. 7. Comparison between Green's function for different a 
and times given in the legend; r is calculated from (|5U|) for 
each a separately. The additional description is the same as 
in Fig. H 



parameters. We are going to find the conditions which 
ensure that the probability densities of FPT and their 
cumulative functions will be the same over the long time 
limit. For t 3> max{tsM,tp} the agreement condition 
reads 



F S M(t;x Ql x M ) = F F (t;x ,x M ) 



(48) 



which is equivalent to 



RsM(t;x ,x M ) = R F {t]x ,x M ) 



From ([35j>, gT]) and (J48J) we get 



V2 



r(l-a/2) y/T(l + a) s/3F=T\z r 



(49) 



The numerical solution to Eq. (|49| has a good approxi- 
mation in the following form (see Fig. [5]) 



r = 3.008a° - 5.471a 4 - 
-0.869a 2 + 0.101a 



3.768a J 
I- 0.463 . 



(50) 

We should note that function ([517)) gives r — > 1 when 
a — > 1~; this result is to be expected since for r = a = 1 
all models provide Green's function for normal diffu- 
sion. The functions GsM{x,t;xo) and Gp-(a;,i; xq) are 
also similar to each other for r given by Eq. (|50|) . which 
is shown in Figs. [6] and [7] In Fig. [6] Greens' functions 
for the Sharma-Mittal and the fractional models calcu- 
lated for a — 0.7 and various times are presented, in 
Fig. [7] for different a and t. Obviously, this similarity 
and the fact that the functions provide the relations ([2]) 
and (|48[) does not necessarily imply any equivalence be- 
tween Greens' functions, e.g. for fixed t and for x — > ±oo 
GsM(x,t;x ) —> (l/|x|) 2 /( 1_r ', and while the asymptotic 




FIG. 8. Comparison between the Greens' functions obtained 
from different models. 



form GF{x,t;xo) ~ |x|exp(— a^ 1 /^ 01 / 2 ^) (a is a posi- 
tive constant) has an exponential character and depends 
on a 0. In Fig. E we can observe that the Greens' func- 
tions of the Sharma-Mittal and the fractional models are 
similar for r, which are calculated from (1501) . in contrast 
to the Gauss and Tsallis models, which do not fit each 
other. This is rather obvious, since the subdiffusive Tsal- 
lis model corresponds to r > 1, whereas the Gauss one is 
only applicable for r — > 1 ~ . 
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V. FINAL REMARKS 

We have compared two models describing subdiffu- 
sion: the model based on Sharma-Mittal nonadditive en- 
tropy, which provides the nonlinear subdiffusion equation 
with derivatives of a natural order and the one based on 
Continuous Time Random Walk formalism which pro- 
vides the linear differential equation with fractional time 
derivative. Sharma-Mittal formalism contains three pa- 
rameters q, Qsm and r whereas the fractional model is 
determined by two parameters a and D a . 

The main results are as follows: 

1. We have shown that Green's function obtained 
within Sharma-Mittal entropy formalism fulfils the 
relation @ (this relation is obvious within the frac- 
tional model) only if the parameters q and a fulfil 
the relation ([23]) and Qsm, D a , a, r fulfil the rela- 
tion (pM|) . We note that relation (IM|) is ambigious 
if the relation between r and the other parameters 
remains unknown. 

2. We have also shown that the first-passage time 
probability distribution F and the probability R 
that a particle has not reached the arbitrary cho- 
sen point xm until time £, calculated within the 
Sharma-Mittal model for r given by the formula 
(l50l) (the other parameters fulfil the relations (|23"]) 
and (|24p). coincide over the long time limit with 
the ones obtained from the fractional model. This 
coincidence is independent of the subdiffusion coef- 
ficient D a . The functions F and R obtained from 
the Sharma-Mittal model for r > 1 (which includes 
the Tsallis model) and for the Gauss entropy model 
(which corresponds to the limit r — > 1~ of Sharma- 
Mittal functions) do not coincide with the results 
of the fractional model. 

The conditions [1] and [5] have a more general char- 
acter since a lot of characterisitcs of the system de- 
rived within both models will also coincide with each 
other, as, for example the time evolution of an amount 
of substance leaving a semi-infinite medium occupy- 
ing the region (— oc,Xm) (then the surface located at 
xm can be treated as a fully absorbing membrane) 
M(t;x M ) = /^C(x o ,0)[l - R(t;x ,x M )]dx , or the 

flux J(t;x M ) = J^C(x o ,0)F(t;x o ,x M )dxo flowing 
through a thin membrane located at xm- Let us note that 
under conditions [T] and [5] the 'similarity' of the Greens' 
functions also occur. This 'similarity' is not here a strict 
definition. It means that the plots of these functions are 
so similar that we cannot decide in an unambigiuos way 
which of these fits better the experimental data. As we 
mentioned earlier, Green's function can be interpreted 
as a normalized concentration of a large number of par- 
ticles starting from the same point at the initial time, 
thus we can suppose that this 'similarity' will occur for 
concentrations with various initial distributions. 



When r is not given by the formula (|50|) (which in- 
cludes the cases r > 1 and 1/3 < r < 0.463), the Greens' 
functions derived within nonadditive entropy formalism 
fulfil the relation © and Eqs. {23), {Ml), © arc 

still valid, but r appears to be independent of a and 
D a . As we have shown in (2(|, subdiffusion described by 
the Sharma-Mittal entropy formalism has a stochastic 
interpretation for 1/3 < r < 1. Namely, it is a pro- 
cess described by the generalized Langevin equation, in 
which the strength of the random force is disturbed by 
the external noise described by the Gamma distribution; 
then, the parameter r is controlled by the mean-value 
of the Gamma distribution. When the parameters of 
both models are related by Eqs. {23]), {MJ) and ([5D|. the 
internal and external noise provides 'similar' effects as 
the noise generating random walk within the fractional 
model. When 1/3 < r < 0.463, the random walk gives 
the effect that the random walk is more hindered by the 
external fluctuations than the random walk described 
by fractional model (a and D a are the same for both 
models). For r > 1, the support of Green's function 
(|2~T1) is finite. In consequence, the subdiffusion process is 
achieved if the velocity of the borders decreases over time 
according to formula {28]) . Let us note that there are two 
mechanisms of the subdiffusion process. The first one de- 
pends on the motion of the borders which limit the area 
penetrated by the walker, the second is the process of 
the walker's movement within the allowed region. Since 
Greens' functions with finite support and nonzero vari- 
ance cannot describe an infinitely divisable process [52j . 
the second mechanism cannot be considered as a random 
walk process in contrast with the fractional model. For 
r — > oo Green's function approaches a constant distribu- 
tion |53j . thus r can be interpreted as the 'measure' of 
the deviation of the probability density of constant dis- 
tribution. Therefore, r represents the 'uncertainly' of 
finding the walker in the near-border regions. 

A potential application of the Sharma-Mittal model 
for r > 1, is an animal searching for food. A simi- 
lar problem was considered in 44|, where the diffusion 
searching process was assumed to be at two stages. The 
first stage consists of searching for a diffusive type near 
the point where the animal stopped its 'quick movement' 
in order to find food, the second stage is the relatively 
fast movement made by the animal in order to change its 
searching region. This movement is assumed to be bal- 
listic, e.g. in constant motion. Changing these assump- 
tions, we get the model which — at least in our opinion 
— can be described by the Sharma-Mittal Green's func- 
tion (|27[) . Namely, we assume that the animal movement 
in the restricted region (stage 2) is carried out with a 
velocity decreasing over time according to formula (|28l) . 
This diminishing velocity is connected with the animal's 
energy, which is gradually lost, especially in the system 
of a comlpex structure. In stage 1, searching inside the 
restricted region can be governed by some specific mech- 
anisms, which are not, in fact, random walk. 

Now we shall turn our attention to discussing a method 
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of extracting the subdiffusion parameters occuring in a 
nonadditivc entropy model from experimental data (a 
similar problem was considered in [54J). The parameters 
q and Qi can be calculated from Eqs. ([23]) . (IM1) . (j2"5j) and 
(f26|) . The parameters a and D a are defined by Eq. ((2J, 
but this equation does not lend itself to experimental ass- 
ignment. Therefore, one needs to use other functions, 
which can be measured experimentally. One of such func- 
tion is M(t;XM) and J(t;XM) as defined above. When 
r is given by ([50]). we get M SM (t;x M ) = M F (t;x M ) and 
Jsm^Xm) = Jf^xm) over a long time limit. Thus, 
assuming that 0,463 < r < 1, each measurement of M 
and J gives the parameter values which can be taken 
from the fractional model. The fractional model gives 
the parameters using a different method, for example, 
the time evolution of a near- membrane layer [l6l [ItJ ■ To 
illustrate this we will now calculate the parameters oc- 
curing in Sharma-Mittal entropy for the subdiffusion of 
glucose and sucrose in gel (1.5 per cent water solution of 
agerose). The subdiffusion parameters found in [ltl [l7| 
are: a = 0.90±0.01, D .90 = (9.8±1.0) x 10" 4 rnm 2 /s° - 90 
for glucose, and A).90 — (6.3 ± 0.9) x 10~ 4 mm 2 / s a9 ° 
for sucrose. Taking r = 0.78 calculated from ([50)1 
and using Eqs. ([23]). (l24l and (j44j) we get q = 1.22 
and Qsm — 6.43 x 10~ 4 mm 2 - 22 /s for glucose, Qsm — 
3.93 x 10 _4 mm 2 ' 22 /s for sucrose (we have omittted here 
the error calculations). 

To extract the subdiffusion parameters from experi- 
mental data for the Sharma-Mittal model when r > 1, 
it is possible to measure the velocity of the borders v(t), 
which allows one to determine the parameters a and B. 
We remark that a similar approach was used to find the 
parameters of the cut-off distribution of postural sway; 
the estimation of the border locations allows us to cal- 
culate the distribution parameters (53[. However, in this 
case, the process was not considered as a stochastic one 



continuously changing over time, but only two arbitrary 
chosen 'times' (postural sway for the old and young) were 
taken into account. 

Summarizing the above remarks, the obvious conclu- 
sion results from the considerations presented in our pa- 
per: subdiffusion cannot be fully characterized by rela- 
tion alone. This remark partially agrees with the 
one presented in (l5| . where the authors concluded that 
relation ([2]) is not enough in order to designate the pro- 
cess as subdiffusion, but its appropriate stochastic inter- 
pretation is also needed. However, taking into account 
the fact that there are situations in which the stochastic 
model has not yet been found, but the models without 
the stochastic interpretation (such as the ones based on 
nonadditivc entropy) are applicable in describing a pro- 
cess, we assume that both the stochastic as well as the 
non-stochastic models can be used to describe subdiffu- 
sion. Under this assumption it seems to be reasonable to 
take relation @ as the definition of subdiffusion under 
the conditions that this relation does not provide another 
important characteristic of the system and its stochastic 
interpretation may not be obvious. In this context the 
Sharma-Mittal model appears to be a relatively universal 
model of subdiffusion where its interpretation and impor- 
tant characteristics (such as first-passage time distribu- 
tion) depends on the parameter r. For special chosen r 
this model gives very similar results to the ones provided 
by the fractional model. 
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